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Abstract 

The time-frequency content of a signal can be measured by the 
Gabor transform or windowed Fourier transform. This is a function 
defined on phase space that is computed by taking the Fourier trans- 
form of the product of the signal against a translate of a fixed window. 
The problem of finding signals with Gabor transform that are maxi- 
mally concentrated within a given region of phase space is discussed. It 
is well known that such problems give rise to an eigenvalue problem for 
an associated self-adjoint, positive concentration operator that has its 
spectrum contained in the unit interval. In this paper, the asymptotic 
behavior of these eigenvalues as the concentration region gets large is 
studied. The smoothness of the eigenfunctions is also examined. 

1 Introduction 

A time- frequency shift of a square integrable function f(t) G L 2 (R) is de- 
fined by: 

p(r, cr)f(t) = exp(-7rzTcr) exp(27ricri)/(t + r). 

The first factor is not essential and is present in order to simplify many of 
the susequent formulae. Ignoring this, the formula clearly has the structure 
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of a time shift of r units together with a frequency shift of a units. The 
Gabor transform or windowed Fourier transform depends on the choice of a 
window function (j)(t) € L 2 (R) and is defined by 

S<j,f{T,a) = (f,p{T,o)4>) 

= J f(t)p{T,a)4>{t) dt 

for any / € L 2 (R). Here (•, •} denotes the standard inner product on L 2 (R). 
The Fourier transform of f(t) is defined by 

f(a)=Tf(a) = J f(t)e- 2mt °dt. 

One easily checks that 

Fp(T,a)f = p(-a,T)f. (1) 

The following proposition collects some formulae that will be useful in 
the sequel. 

Proposition 1 Let f,(/>E L 2 (R). 

1. JJ\S^f(r,a)\ 2 drda = \\f\\ 2 Ur. 

2. <S^/(r, a) = cr, r) where f and <f> denote the Fourier transform. 

Proof: The first item is a well known special case of Moyal's identity. See 
P]. The second item is proved by using PlanchereFs theorem together with 
equation [l|: 

S<t>f{T,o) = (f,p(T,a)(j)} 

= {Ff,Fp{T,o-)<i>) 

□ 

Let n c R 2 be a region within the time-frequency plane. For the re- 
mainder of the exposition, will be assumed to be normalized to have L 2 
norm equal to one: \\4>\\ 2 = 1. We will measure the energy of / contained 
within f2 by 

£u(f) = J l^f(r,a)\ 2 drda. 



2 



Associated to the energy is the Hermitian symmetric form defined by: 

£n(f,g)= / / {f,p{r,cr)(j)){p(T,a)(j),g)dTda 
J Jn 

for any f,g£ L 2 (R). Proposition [l] insures that £ is bounded in the follow- 
ing sense: 

£n(f,9)< ll/ll \\9l (2) 

A standard duality argument produces a bounded, positive, self-adjoint op- 
erator C n : L 2 (R) -> L 2 (R) such that 

£n(f,g) = (Caf,g) 

for all /, g £ L 2 (R). Let i^i,^,- ■ ■ be a complete orthonormal basis of 
L 2 (R) and Ai, A2, • • • be a decreasing sequence of nonnegative real numbers 
for which 

Cni>k = Afc-^fc- 

The bound in equation |2| forces 1 > Ai > A2 > • • •• 

It is clear that the unit norm signal with the maximum energy within f2 
is ipi. Moreover, this maximum energy is given by Ai. More generally, the 
quantity 

n(A) = card{A; : A& > A} 

determines the maximum dimension of a subspace V E L 2 (R) of signals for 
which 

^(/)>A||/|| 2 

for all / G V. Thus the distribution of these eigenvalues as well as the 
smoothness of the corresponding eigenvalues are of interest within the con- 
text of signal processing. 

2 A Reproducing Kernel Hilbert Space 

It will be convenient to analyse the distribution of the eigenvalues using the 
reproducing kernel Hilbert space setting described in this section. 

Proposition [l] implies that Sa, is an isometry of L 2 (R) onto a closed 
subspace C L 2 (R 2 ). Let be the orthogonal projection from L 2 (R 2 ) 
to Vcj). Let F G and / = S^ l F. Then, applying the Cauchy-Schwarz 
inequality, one has 

l^(^-o,^o)| = (f,p(T ,a )<P) < ll/ll ||p(7b,<roMI = ll^ll- 
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Therefore point evaluation at a given (to, ctq) of functions in V<j> is a bounded 
linear functional. Let Kr TOjao \(r, a) = K(To,ao]T,a) denote that element of 
V<p for which 

F(r Q ,a ) = J J F(r,a)K 

In order to identify the kernel explicitly one only has to notice that 

F(r ,cr ) = (S^F^faaoW). 
Since is an isometry, 

F(r ,cr ) = {F,S<j,p{To,oo)4>) 

= J J F(T,a){p(T ,a )(j),p(T,a)(f>) drda 

Consequently, 

K(t ,(To;t,(t) = (p{T ,a )4>,p{T,a)(p). 
Finally, if F G Vi we have 



J J F(t, a)K (nuao) (T, a) drda = 



since K/ T0(T0 \(t, a) G V^. These remarks are summarized in the following 
proposition. 

Proposition 2 For any F G L 2 (R 2 ) the projection V^F is computed by 
P^F(t,<t) =J J F(r',a')K {T!a) (T',a')dr'da'. 

Our next task is to use the isometry to identify the concentration 
operator Cq explicitly. In particular, the next proposition provides a formula 
for S^CqS^ 1 on V^. 

Proposition 3 For any F G V<j>, 

S^CnS^F^a) = J J F{t> \a')K {T ^{r' ,a')dr'da' '. 
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Proof: Let S^f = F and S^g = G. Then 

£<n(f,g) = {f,p(r,cr)<j)){p(a,T)(j),g)dTda 



J J^S^fS^g drda 



FG drda 

<> 

F\ n GdTda 

= J J V4F\ n )Gdrda. 

On the other hand, £n(f,g) = (Cq/,<?). Putting this together easily gives 
the formula in the proposition. □ 

In view of proposition ||, the spectral properties of the concentration 
operator are identical to those of the operator from to itself defined 
by 

where xn is the operator that restricts functions to the set 0. 

Since the function (j) is fixed in the ensuing discussion, its use as a sub- 
script will often be supressed in the sequel. 



3 Distribution of the Eigenvalues 

In this section we examine the behavior of the eigenvalues as the concen- 
tration region 0, gets large. In particular, we show that the number of 
eigenvalues that cluster around one grows like the area of the concentra- 
tion region. Both the philosophy and proofs in this section bear a great 
resemblance to the work of H.J. Landau in j?], ||]. 

Theorem 1 Let $7 be a bounded measurable set and 

n(A, f2 r ) = card{/c : Afc(S7 r ) > A} 

where f2 r = {(rr, ra) : (r, a) £ Q}. Then 

hminf^4 = l 
r— >oo area(ii r ) 

as r — > oo. 
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The essential idea behind this theorem, as in the bandlimited case treated 
by Landau, is that J2 an d J2 \^k\ 2 grow at the same asymptotic rate as 
r — > co. The following lemma enables us to work with the operator xvlF on 
L 2 (Q) as opposed to the operator Vxn on V. 

Lemma 1 The operator xtiP on L 2 (Q,) and the operator Vxo. have the 
same nonzero eigenvalues with multiplicity. 

Proof: Let F E V be a nontrivial eigenfunction with a nonzero eigenvalue of 
the operator Vxn'- FXnF = A-F. Since both F and A are nonzero, XttF ^ 0. 
On the other hand, applying the restriction operator to both sides we have 
XciFxnF = XxnF. A similar argument shows that every nonzero eigenvalue 
of XnP is a lso an eigenvalue of Vxn- 1=1 



Proof: By lemma |], it suffices to consider the trace operator {xn r F) 2 - The 
kernel of this operator on L 2 (fi r ) is determined by the formula 



(xn r V) 2 F(r,a) = 

/ / / In r xn r Hr", a")K {T , y) (r" , a")K^ a) (r', a') air" da" dr' da' . 



The kernel is continuous on the bounded set VL r x Q, r so by Mercer's theorem 
the trace of (xvirF) 2 is 



Lemma 2 




1 as r — > co. 



//„ 




K(T,a;T',a')\ 2 dTdadr'da'. 



First note that 



K{r,a-r',a') 



(p(r, a)(f),p(T',a')(fi) 

(p(- t ',-v')p(i-, <t)4>,4>) 

exp (7n(rcr — ar )) (p(r — r , a 

exp (m^ra' — ar')) H{t — r', cr 




Therefore, it is enough to estimate 



// If \H(T-T',a-a')\ 2 dTdadT'da'. 



By proposition [T] the 1? norm of H is one. The theorem then follows from 
the following lemma. □ 
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Lemma 3 Let f G L 1 (R n ) wiift /(x) > and J f(x)dx = 1. Let QcR" 
is a bounded set of positive measure. Then 



n [ I f( x -y)dxdy^\Q\ 

JrQ JrQ 



as r — > oo. 



Proof: Apply the following change of variables: 

x = £ + rrj y = rrj 

to get 

r~ n [ f f(x-y)dxdy= [ f f(0^d V . (3) 

JrQ JrQ Jq Jr(Q-t]) 

Recall that almost every point of a set of positive measure is a point of 
density: 

lim \ B ^^®\ = i for almost all r, G Q. 
Hence for almost all rj G Q and arbitrary R > 0, 

HQ- v )nB R (o)\ 

ttR 2 ^ 
as r — > oo. It is then straight-forward to argue that 



/ /(Ode - / /(Ode 



/r(Q-r?)nB fl (0) ^Bjj(O) 

as v — ► oo for almost all rj e Q. Since 

/ /(Ode < / /(Ode < i 

it follows that, for almost all ?? G Q 

( /(Ode < hminf f < lim sup / /(O^ < 1- 

JS fl (0) r ^°° Jr{Q-n) r^oo Jr(Q-r)) 

Since i? > was chosen arbitrarily, it follows that 



/ /(Ode = i 

Jr(Q-n) 



lim 

r ^°° -/r(Q-7/) 
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for almost all rj E Q. The result then follows from Lebesgue's dominated 
convergence theorem. □ 

PROOF of Theorem |l|: By Mercer's theorem, the trace of the operator 
T>X is 

Afc = II K(t, a, t, a) drda 

(p(r, a)<j),p{T,o)4>) drda 
= area($7 r ). 

Choose r so large that |Afc(O r )| 2 > (1 — e)area(f2 r ). Then 
(1 - e)area(O r ) < ^\\ k (n r )\ 2 

n 

k=l k>n 
n 

< (1 — A) Afc + Aarea(f2 r ) 

k=l 



As a consequence, 



< (1 — X)n + Aarea(O r ) 



1 — A — e 

n > areaiiL). 

~ 1-A y J 



liminf > 1. 



Since e > was chosen arbitrarily, we have that 

ra(A, f2 r ) 
area(ri r ) 
□ 

The following refinement of lemma ||, in the case when Q is a bounded 
domain with C 1 boundary and / has sufficiently strong decay, is useful in 
order to estimate the size of the plunge region. 

Lemma 4 Let Q be a domain with C 1 boundary and f(x) as in lemma |^. 
Moreover, assume that 



C 

j{x) dx\ < 

for some constants p, C > 0. Then the error term 



1-/ f{x)dx\<- (4) 
J BAG) rP 



\l I f(x-y)dxdy-\Q\ 

r n JrQ JrQ 
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is 0(^) as r —> oo. 

Proof: Given a 5 > 0, let Qs = {x G Q : dist(x, dQ) > 5}. The second 
iterated integral in equation |3| can then be expressed as follows: 

5 ° [ I r ( V )dH( V )dp (5) 

JdQp 
where 



Jr(Q-rj) 

dTL is Hausdorff n — 1 dimensional measure, and 5q is any positive number 
larger than the diameter of Q. Using equation |] to estimate I r {f]) yields 



f° [ I r ( V )dH( V )dp-\Q 1/r \ < [ 5 ° [ \I r ( V )-l\dH( V )dp 

Jl/r JdQp Jl/r JdQp 



So 

fr JdQ p 



< C f° [ J—dH^dp. (6) 
Jl/r JdQp {rp) p 

Using the fact that the Hausdorff n — 1 measure of the sets dQ p are uni- 
formly bounded, one easily has an estimate of the required type for the first 
expression in inequality ||. Moreover, since < I r (ri) < 1, 



Io /r JdQp J r(v) dH(rj) d P ~ meas(Q\Q 1/ , 

-1/r 


/o 1A 7a Qp dH( V )dp 



Io /r IdQp (l-I r (v))dH( V )dp< (7j 



< meas(Q\Qi/ r ) < C/r 

Putting equations |6| and |7| together via the triangle inequality yields the 
sought after result. □ 

Theorem 2 Let Q, be a bounded domain with C 1 boundary and 

n(\,p,£l r ) = card{/c : Afc(S7 r ) G [A,//]} 

where < A < p < 1. Moreover, assume that the function \S<p(j){T, o~)\ 2 
satisfies condition 0. Then, there is a constant C > for which 

n(A, p, Q r ) < Cr. 
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Proof: Mercer's theorem applied to the positive operator ~ (xn^) 2 

yields that 



^Afc(Q r ) -X k (n r ) 2 = |Jlr| -J J J J \H(t -t' \<j -a')\ 2 drdT dodo '. 

Applying lemma Q yields 

0<^A fc (^ r )-A fc (fi r ) 2 <Cr. 
Therefore, if e G (0, 1/2) then 

n(e, 1 - e, Q r )(e - e 2 ) < ^ A fe (O r ) - A fc (ft r ) 2 < Cr. 
The result follows. □ 

4 Regularity of the Eigenfunctions 

In this section we will be concerned with the regularity of the eigenfunctions 
of the concentration operator Cq. In order to do this it will be necessary to 
work with the integral kernel of this operator. 

Proposition 4 The kernel kn(x , y) of the operator Cq is 



kn(x,y)= / / p(T,o)(p(x)p(T,o)4>(y) drdo. 



n 



Proof: Using the definition of the concentration operator, one has 

(Cnf,g) = £n(f,g) 

(f, P(r, <y)4>) (p{r, o)4>, g) drdo 

n 

f(y)[ / p(r,o)(j)(x)p(T,o)(p{y) drdo) g(x)dxdy. 



The result follows. □ 

Suppose that 7(5) is a positive decreasing function defined on [0, 00) 
satisfying the following conditions: 

• 7(0) = and lim^oo 7(5) = 0, 
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• there is an eo > such that / \^(s)\ p ds < oo for all p € (2 — eo,2] , 
and 

• for any sq > 0, j(s — so) = 0{ r ){s) 1_e ) for all sufficiently small e > 0. 
Examples of functions j(s) that satisfy such conditions include 

• 7( s ) = (1 + s 2 )~ q / 2 where q € (1, oo), and 

• 7(s) = exp(— n\s\ q ) where n, q > 0. 

Proposition 5 If \4>(t)\ < j(\x\) for all t then any eigenfunction with 
nonzero eigenvalue ip of the concentration operator Cci satisfies an estimate 
of the form 

mt)\ < cm) 1 -* 

for all t. 

Proof: Suppose that Cnip = Atf) with A > 0. The first step is to estimate 
the kernel, kfi(x,y): 



\kn(x,y)\ 



e 2na( - x - v U(x + T)<j){y + t) drda 

n 



< C 7 (||x|-fl[)7(||!/|-fl|) 

< CTdx^-Tdyl) 1 "* 

where R = sup{|(r, er)| : (r, cr) € Vi}. Consequently, using the Cauchy- 
Schwarz inequality, we have the estimate: 

A|V(x)| < CiW) 1 -' j -yiW-'^dy 
< Cj(\x\) l - C . 

□ 

Let Cl = {(t, a) : (— a, r) £ The second item in proposition [j] implies 
that: 



£a(f,g) = J J^s^fS^g 



n 



In particular, 

and V'l) ' " " are tire eigenfunction of the operator Cq. Putting this together 
with the preceding proposition yields the following theorem. 

Theorem 3 If 

\<t>(t)\ < 7l(M) and 

< 72(H) 

where both 71 and 72 satisfy the conditions described above, then 

\^(t)\ < Cjxdtl) 1 -" and 
\fa)\ < C l2 (\a\) l ~\ 

Various special cases merit mention. 

• Let 

Mi\ _ / 1 " 1*1 if 1*1 < 1 1 

ni> ~\ otherwise. J 

Then kn(x,y) = whenever \x — y\ > 2. Hence, any eigenfunction, 
with nonzero eigenvalue, must have support contained in the interval 
[—3, 3]. Moreover, the Fourier transform of cj> is the Fejer kernel which 
is bounded by 7(5) = (1 + s 2 )" 1 . As a consequence, if Hence, if ip is 
any eigenfunction with nonzero eigenvalue, 

support(^) C [—3, 3] 

fa) = 

In particular, ijj' is square integrable. 

• If is Schwartz class, one can verify directly using proposition |4| that 
the kernel &q is also Schwartz class. A straightforward interchange of 
the differentiation and integration symbol then shows that tp is also 
Schwartz class: 

\\t a Dfy(t)\ = \ft a D^k n (t,s^(s)ds\ 

< \\t a D?k n (t,s)\\ 

Since the kernel is Schwartz class, the last expression is bounded inde- 
pendently of t. Thus, if 4> is Schwartz class, so are any eigenfunctions 
associated to nonzero eigenvalues. 
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• If 4>(t) = exp(— Ct 2 ), theorem || implies that the eigenfunctions of 
the concentration operator with nonzero eigenvalue are C(exp(— C(l — 
e)t 2 )). Moreover, since 0(a) = y/ir/C exp(— ir/Ca 2 ), the Fourier trans- 
forms of these eigenfunctions are C(exp(— ^(1 — e)cr 2 ). In particular, 
these functions are entire by the Paley- Wiener theorem §. Further- 
more, by differentiating under the integral sign (as above, and in the 
formula for the kernel in proposition |4j), one can show that all deriva- 
tives of these eigenfunctions satisfy the same decay estimates. 

Theorem |3| is useful in that it provides quantitative bounds on the decay 
(and eventually regularity) of the eigenfunctions for a wide class of eft. Among 
the works in the literature that bear on the problems discussed here, we 
remark that Pietsch Q provides estimates on the decay of the eigenvalues 
for given decay and regularity of the integral kernel. In particular, his results 
imply that when 4> is in Schwartz space, the sequence of eigenvalues decays 
faster than the inverse of any polynomial in the index n. Janssen § obtains 
the same estimate in the Schwartz class case. 

For <p(t) = exp(— Ct 2 ), the work of Daubechies Q is relevant. She solves 
the eigenvalue problem explicitly for circular (or elliptical) domains, and 
obtains the Hermite functions as eigenfunctions. She also obtains a formula 
for the eigenvalues, which decay exponentially in the index n. Our approach 
shows that for any bounded measurable domain, the eigenfunctions are ana- 
lytic and have quadratic exponential decay on the real line. Since the kernel 
is Schwartz class, || implies that the eigenvalues decay faster than the in- 
verse of any polynomial in the index n. In fact, a stronger statement can be 
obtained using a result of Janssen 0] . Note that the integral kernel is 

kn(x,y) = j J n exp(2iria(x — y))exp(—C((x + r) 2 + (y + t) 2 )) drda. 

By expanding the exponentials in power series and integrating term by 
term, it is easy to show that this kernel is analytic in each variable, and 
satisfies the bound 0(exjp(-ARe[x] 2 + BIm[x] 2 - CRe[y] 2 + DIm[y] 2 )), for 
some A, B,C,D > 0. Since the kernel is positive definite as well, it follows 
from theorem A.l in || that the eigenfunctions belonging to nonzero eigen- 
value are analytic and satisfy 0(ex.p(—ARe{x] 2 +BIm{x] 2 )), some A, B > 0. 
At the same time, the eigenvalues are in C(exp(— no)), some a > 0. 
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